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We study the dielectric function of the homogeneous semiconductor hole liquid of p-doped bulk 
III-V zinc-blende semiconductors within random phase approximation. The single-particle physics 
of the hole system is modeled by Luttinger's four-band Hamiltonian in its spherical approxima- 
tion. Regarding the Coulomb-interacting hole liquid, the full dependence of the zero-temperature 
dielectric function on wave vector and frequency is explored. The imaginary part of the dielectric 
function is analytically obtained in terms of complicated but fully elementary expressions, while 
in the result for the real part nonelementary one-dimensional integrations remain to be performed. 
The correctness of these two independent calculations is checked via Kramers-Kronig relations. 

The mass difference between heavy and light holes, along with variations in the background 
dielectric constant, leads to dramatic alternations in the plasmon excitation pattern, and generically, 
two plasmon branches can be identified. These findings are the result of the evaluation of the full 
dielectric function and are not accessible via a high-frequency expansion. In the static limit a beating 
of Friedel oscillations between the Fermi wave numbers of heavy and light holes occurs. 



PACS numbers: 71.10.- 



71.10.Ca, 71.45.Gm 



I. INTRODUCTION 

The interacting electron gas, combined with a homo- 
geneous neutralizing background, is one of the paradig- 
matic systems of many-body physicsi"^. Albeit the result 
of drastic approximations, its predictions provide a good 
description of important properties of three-dimensional 
bulk metals and, in the regime of lower carrier densities, 
n-doped semiconductors where the electrons reside in the 
s-type conduction band. 

On the other hand, in a p-doped zinc-blende III-V 
semiconductor such as GaAs, the defect electrons or holes 
occupy the p-type valence band whose more complex 
band structure can be expected to significantly modify 
the electronic properties. Moreover, the most intensively 
studied ferromagnetic semiconductors such as Mn-doped 
GaAs are in fact p-doped with the holes playing the key 
role in the occurrence of carrier-mediated ferromagnetism 
among the localized Mn magnetic moments^. Thus, such 
p-doped bulk semiconductor systems lie at the very heart 
of the still growing field of spintronics^, and therefore it 
appears highly desirable to gain a deeper understanding 
of their many-body physics. 

Ab-initio-type approaches to the description of ferro- 
magnetic semiconductors constitute an important sub- 
field of this endeavor, and there is a lively discussion on 
strengths and weaknesses of the pertaining concepts and 
numerical techniques^. In the present paper we will fol- 
low a somewhat different route by developing an analyt- 
ical theory of the most prominent class of host materials 
given by p-doped bulk III-V zinc-blende systems such 
as GaAs. Specifically, we investigate the dielectric func- 
tion of the interacting hole liquid within random phase 
approximation (RPA)2i2] where the non-interacting hole 
system in the valence band is described by Luttinger's 
Hamiltonian in the spherical approximatiorP". We evalu- 



ate the zero-temperature dielectric function in the entire 
range of wave vectors and frequencies building upon a re- 
cent study where the problem was analyzed in the static 
limit, and in the case of large frequencies^. Another pre- 
vious work investigated, among other issues, properties 
of Hartree-Fock solutions of the two-component carrier 
system consisting of heavy and light holes^l. Moreover, 
very recently Kyrychenko and Ullrich have put forward 
a study of holes in magnetically doped III-V systemJISl 
by modeling the band structure by an 8 x 8 fc • p Hamil- 
tonian (similar as in the present work) while disorder ef- 
fects and interaction among the carriers are treated by a 
combination of equations-of-motion techn iques and time- 
dependent density functional theorjli^liil. Further below 
we will compare our fairly analytical results with the ones 
of Ref.'^ which rely heavier on numerical evaluations. 

Finally we mention a series of related recent studies 
of the dielectric properties of two-dimensional fermionic 
systems (instead of three-dimensional bulk semiconduc- 
tors) whose single-particle states carry a non-trivial 
spinor structure. These include n-doped quantum wells 
with spin-orbit couplingiSES and two-dimensional hole 
systemalS. Other recent investigations have dealt with 
planar graphene sheets where an effective spin is incor- 
porated by the sublattice degree of freedorrP^Iill 

The plan of the paper is as follows. In section |TT] 
we give an overview on elementary properties of the 
single-particle Hamiltonian describing the band structure 
around the F-point, and on the many-body formalism 
leading to the RPA result for the dielectric function. In 
Section III we present explicit analytical expressions for 



the free polarizability; the corresponding derivations are 
deferred to the appendices. Section |IV] discusses physical 
properties of the dielectric function and its full depen- 
dence on wave vector and frequency. Special attention is 
paid to the static limit, and to the limit of large frequen- 



cies. We close with conclusions and an outlook in section 
lYl 



II. PRELIMINARIES: HAMILTONIAN, 

EIGENSYSTEM, AND MANY-BODY 

FORMALISM 



Luttinger's Hamiltonian describing heavy and light 
hole states around the F-point in III-V zinc-blende semi- 
conductors reads in its spherical approximatiorP, 
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Here mo is the bare electron mass, p is the hole lat- 
tice momentum, and S are spin-3/2-operators, result- 
ing from adding the I = 1 orbital angular momentum 
to the s = 1/2 electron spin. The dimensionless Lut- 
tinger parameters 71 and 72 describe the valence band 
of the specific material with effects of spin-orbit coupling 
being included in 72. We note that , while the present 
work is mostly motivated by III-V semiconductors, the 
above model for the Fg valence band also applies to other 
systems with zinc-blende or diamond structure including 
elemental semiconductors like Si and Ge, but also zero- 
gap semiconductors such as HgSe and HgTe. 

The above Hamiltonian is rotationally invariant and 
commutes with the helicity operator X = (k- S)/k, where 
k = p/h is the hole wave vector. The heavy (light) holes 
correspond to A = ±3/2 (A = ±1/2) with the energy 
dispersions 
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and heavy {h) and light {I) hole masses m^^i = mo/(7i T 
272). The corresponding eigenstates are given by 



and the remaining eigenspinors |X-3/2(fc))7 lx--i/2(fc)) 
can be obtained from the above ones by spatial inversion 
■di-^Tr — 'd,ip'-^ip + Tr. In what follows, mutual overlaps 
squarecP' between spinors will be of key importance, 



l(X|(fci)IX|(^2))|' 
|(xi(fci)|xi(fe))P 
|(X3(fci)|xi(fe))P 



Combining the above single-particle Hamiltonian with 
Coulomb repulsion among holes and a neutralizing back- 
ground, the dielectric function within RPA at wave vector 
q and frequency w is given by ^ - 
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Here V{cj) is the Fourier transform of the interaction po- 
tential, and the free polarizability reads 
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where V is the volume of the system. Using the con- 
ventional basis of eigenstates of 5'^ and introducing po- 
lar coordinates k — fc(cos(psin'i9,sin(^sin'i9,cos^) , the 
eigenspinors \x\{k)) of the helicity operator read explic- 
itljP 
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with f(k, A) being Fermi functions. In what follows 
we will concentrate on the case of zero temperature 
and Coulomb repulsion, V{q) = e? / {er£oq^) where e,- 
is a background dielectric constant taking into account 
screening by deeper bands. 



III. THE FREE POLARIZABILITY 



We now present our analytical results for the real and 
imaginary part of the free polarizability. Details of the 
derivations are can be found in appendices] A l|and A 2 A 



discussion of the physical properties of the corresponding 
dielectric function follows further below in section IIVI 



Defining 




Xhh{q,^) = 


1 

(27r)3 



d-'k - 



1 3 



k-{k + q) 



k<k,. 



2 2 k^{k + q)^ 



i{uj + iO)~ (e,,(fc + ^-£,,(fc)) 



ri(a; + iO)+ £h(/s + 9) -£/j(fc) 



(11) 



Xhiiq^i^) 



1 



d-^/c 



3 3 



{k-{k + q)y 



(2^)'A<fc. (2 2 fc2(fc + g)2 

1 



H{u; + iO)- [ei{k + q)-ehik) 

1 

fi(a; + iO)+ (ei{k + q)-ehik)]_ 



, (12) 



one can formulate the polarizability ( 10 ) as follows 
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where the remaining quantities Xuiqi ^) &n.d xih{q, ^) are 
given by ( 11 ) and ( 12 1 via the replacement h -(r^ I, and kh 



(ki) is the Fermi wave number for heavy (light) holes cor- 
responding to the common Fermi energy e^PS! Introduc- 
ing the obvious decomposition Xap{qi'^) = Rapiq,^) + 
ilap{q,uj) with real functions Rap{q,uj), Iap{q,uj) (a,/? G 
{h,l}), we will now analyze the real and imaginary part 
of the free polarizability of the hole gas. The respective 
expressions to be presented below are the result of inde- 
pendent calculations and perfectly fulfil Kramers-Kronig 
relationJiHsl, 



A. The real part of the free polarizability 



Following the steps detailed in appendix |A 1| the real 
part of the free polarizability can be obtained as 
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where (w i— >■ — w) denotes terms with the sign of the fre- and the remaining contribution Rii[q,uj) + Rih{q^Lo) fol- 
quency changed compared to the preceding expression, lows via h ^ I. In the limit ?n/j — mi the first two 



lines in Eq. ( 14 ) express the result for the standard text- 
book case of a fermion gas without spin-orbit coupling^™ 
while all other terms vanish in this limit and represent 
corrections arising from m^ 7^ rn;. The contribution in 
the third line of Eq. ( |l4| ) is constant, i.e. independent of 
(f and uj. However, in the limit of large frequencies this 
term cancels against the terms in the last two lines of 
the above equation such that limi^_>.oo xoi'li'-^) — 0. The 



integral occurring in the fourth line of Eq. ( 14 ) is ele- 
mentary but lengthy (cf. appendix |A 1[) while all other 



integrals cannot be cast be cast into elementary expres- 
sions. Note that in the fifth line of the above expres- 
sion a proper Cauchy principal value (denoted by V) 
occurs. This mathematical detail arises from the Dirac 
identity, and the corresponding integral does for negative 



frequency a; < not converge in the general sense. The 
occurrence of such nontrivial principal values is also a 
technical difference to the standard jellium model. 



B. The imaginary part of the free polarizability 

As the free polarizability xo ('^j t) is a real quantity, 
let us concentrate on non-negative frequencies a; > 0. 
The regions of non-zero contribution Ifiii{q,uj) in the q- 
w-plane are given explicitly in table IT] and are depicted 
for typical system parameters in Fig.na). In region I and 
n Ihh{q,oj) is given by 
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respectively, and are zero for all other values of q and 
UJ. The region boundaries given in table [l] are completely 
analogous to the ones found for a standard jellium gas 
of spinless particles with mass rtih and Fermi momentum 
kh ; for more details see appendix |A2[ The contribu- 
tions to the imaginary part occurring in these regions are, 
however, clearly different from the standard case. The 
regions of nonvanishing contributions of and the corre- 
sponding expressions for In {q, u) can be obtained directly 
via the replacement h 1-^ I. 



The cases of the remaining expressions Ihiiq^u) and 
Iih{q^Lj) are substantially more complicated. It is useful 
to distinguish two separate terms. 
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and likewise for Iih{q,itj). The corresponding regions of 
nonzero contribution to Ii^i{q,uj) and I[^{q,uj) are given 



'/i./v 



in tables ITT] and III respectively, and plotted in Figs, lit) 
and d) for typical parameters. Now defining 
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If^i{q,uj) in regions I and II can be expressed as 

I- ^hiiq^^) = G+{q,uj;kh,kl;mh,mi) , (19) 
II- ^hii^'^) = G+iq,^;kh,kh:'n^h,mi) , (20) 
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TABLE I: Boundaries of regions of nonzero imaginary con- 
tribution Ihh{q,ijj)- The boundaries for Iu{q,uj) are obtained 
via the replacement ht-^ I. 



The nonzero contributions to /^^((f, oj) in regions III and 
IV of table |ll] are given by 

III- ^hii^^^) = G-{q,uj;kh,k'^;mh,mi) , (23) 
IV: Ihiiq,^) = G-{q,uj;kf^,k];;mh,mi) , (24) 

The nonvanishing contributions to Iii^{q,uj) can be ex- 
pressed in a similar manner. For Ii^{q,uj) in regions I 
and II of table IIIII one finds 
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Likewise, the contributions to In^{q,uj) in regions I and 
II are given by 

III: Iih^iq,uj) = G-{q,uj;ki,k'[;mi,mh) , (29) 
IV: Iih^iq,uj) = G-{q,uj;ki ,k^;mi,mh) , (30) 



IV. THE DIELECTRIC FUNCTION 
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TABLE II: Boundaries of regions of nonzero contributions to 
I^j{q,u}) (regions I, II) and I^i{q,Lj) (regions III, IV). 
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TABLE III: Boundaries of regions of nonzero contributions to 
^uM^) (regions I, II) and /,^(g,w) (regions III, IV). 



Let us now analyze the RPA dielectric function result- 
ing from the above free polarizability. We first concen- 
trate on the effect of the mass difference between heavy 
and light holes. To this end we eliminate effects of the 
dielectric background by putting Sr — ^, and we fix 
the total density n — rih + ni, n^/i = fc^;;/37r^, to 
n = O.Olnm^"^. Figs. [2] [3] show the reaP^Jand imaginary 
part of the dielectric function as a function of wave num- 
ber and frequency in a color-coded density plot, whereas 
in the modulus of £^^^{q,u]) is shown. The top left 
panel in each figure illustrates the textbook caseH^: of 
equal masses rriH — mi = mo/2 with its well-known plas- 
mon dispersion U!{q) determined by e^^ {q,Ld{q)) = 0. 
With increasing mass difference between heavy and light 
holes a more complex structure arises and the plasmon 



dispersion splits into two branches as seen in the bottom 
panels of Fig. |4} A branch with comparatively high ener- 
gies at small wave numbers is accompanied by a branch 
at lower energies and large wave vectors. It is an inter- 
esting speculation whether one can interpret these two 
plasmon branches in analogy to phonons: On one branch 
both heavy and light holes possibly perform (speaking in 
classical terms) joint collective oscillations of charge den- 
sity (analogous to acoustic phonons), while on the other 
branch they oscillate opposite to each other (similar to 
optical phonons). We leave this particular issue to future 
investigations. 



Finally, Figs. [5j[6j show the free polarizability as a func- 
tion of frequency at different wave vectors for the same 
choice of heavy and light hole masses as in the previous 




FIG. 



Regions of nonvanishing contributions to a) 



Ihh[q,i^), b): Iu{q,u)), c): I+{q,iu) (solid lines) and 4,(if,Lj) 
(dashed lines), and d): I ^{q ,uj) (solid lines) and Ij^{q,uj) 
(dashed lines); cf. tables [l |lll[ We have chosen the mass pa- 
rameters of GaAs, mh = 0.5mo, mi = O.OSmo, and a hole 
density of n = O.Olnm"^. 
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FIG. 2: The modulus |Re(e^^^(g,a;))| of the real part of the 
RPA dielectric function as a function of wave number q and 
energy fio; for a model system with £,. = 1. The ratio of heavy 
and light mass is varied at constant mu + mi — mo, and the 
total hole density is n = O.Olnm"'^. 



figures. 



Let us now discuss our results for the dielectric func- 
tion with respect to concrete III-V semiconductors. In 
order to make contact to typical ferromagnetic semicon- 
ductor systems^, and to compare with results of Ref.l^ 
we choose here a higher carrier density of n = 0.35nm~'^. 
We consider four typical III-V sys tems whose relevant 
parameterfPil are given in table IV Note that now also 
the background dielectric constant Sr plays a nontrivial 
role. In Figs.[7J|8J[9]we have plotted the reaP^land imagi- 
nary part, and the modulus, respectively, of the dielectric 
function as a function of wave number and frequency. 
As seen from Fig. |9J the zeros of the dielectric func- 
tion e^^ {q,uj) defining the plasmon excitations form a 
clearly more complex pattern than in the standard jel- 
lium liquid, and, as in the previous case, two dispersion 
branches can be identified. In particular, the plasmon 
excitation in GaAs at small wave vector occurs slightly 
below 0.3eV which agrees very well with Fig. 4 of Ref.^^ 
where a more complex model for the band structure was 
used. However, differently from the findings there, we 
can identify two plasmon dispersion branches with small 
damping. Moreover, In Fig. [TO] we show the free polar- 
izability xoC?,'^) as a function of frequency at different 
wave vectors for the same semiconductor systems. Again, 
the imaginary part for GaAs agrees nicely with data given 
in Fig. 2 of Ref.'^i'. In this regime the imaginary part 
of the free polarizability is dominated by transitions be- 
tween heavy-hole states, i.e. the main contributions is 
Ihhiq,!^), in accordance with Ref.'i^'. 
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FIG. 3: The imaginary part lm{e'^^ "^ {q, to)) of the the RPA 
dielectric function as a function of wave number q and energy 
hid for a model system with e^ = 1. The ratio of heavy and 
light mass is varied at constant mjj + mi = toq, and the total 
hole density is n = O.Olnm"^. 



0.7 

0.6 

> 0.5 

S) „ „ 

!n 0.3 



iiiL=0.5niO 


■^' 


-2 


1 
0.7 


^ d 


^ 


1.5 

1 


0.6 
0.5 
0.4 
0.3 


^^P 




0.5 


0.2 


^^^ 







0.1 



0.5 



0.6 1 1.5 2 

wave number [1/nm] 




0.7|^^^^^^^H 
I mL=0.05mO ^^H 



1.5 



0.5 1 1.5 2 

wave number [1/nm] 



■> 

"5-0-1 

¥ 

" -0.2 



-0.4. 




n=0.01nm " q=0.4iim 

m^=0-5mg, m,=0-5mj 
m^=0.65inQ. m|=0.35tn^ 
m^=0.1^mg, m|=0.2mj 
m|_=a95m„, m|=a05m, 




q=0.8iim 




0.1 0.2 

hwl(l%) [eV] 



0.1 0.2 0.3 

h(al(2lt) [eV] 



FIG. 6; The imaginary part of the free polarizability xo(?, ^) 
as a function of frequency at different wave vectors for the 
same choice of heavy and light hole masses as in the previous 
figures. 



FIG. 4: The modulus \e^^'^{q,uj)\ of the RPA dielectric func- 
tion as a function of wave number q and energy hui for a model 
system with e^ = 1. The ratio of heavy and light mass is var- 
ied at constant mn +mi = mo, and the total hole density is 
n — O.Olnm"^. The dark areas indicate zeros of the dielectric 
function corresponding to plasmon excitations. 



> -0.5 




-0.5 



n=0.01nm q=0.4iim 



m|^=0-5mj|, m|=0-5mj| 
m^=0-65ing, m,=0-35inj 
m^=0-8mg, m,=0-2mg 
lu =0-95m„, m,=0-05m„ T 



q=0.6iim 




q=1.0nm 



0.1 0.2 

hco/{2jt) [eV] 



0.1 0.2 0.3 

h(aiait) [eV] 



FIG. 5: The real part of the free polarizability xo('f)'^) as a 
function of frequency at different wave vectors for the same 
choice of heavy and light hole masses as in the previous fig- 
ures. 



Static limit 



In the static limit cj = 0, the rather complex contri- 
butions ( 14 1 to the free polarizability of the hole system 
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FIG. 7: The modulus |Re(e^^^((f,a;))| of the real part of the 
RPA dielectric function as a function of wave number q and 
energy huj for various semiconductor systems at a total hole 
density of n = 0.35nm~^. 



AlAs 
GaAs 
InAs 
InSb 



mo 



0.47 
0.5 
0.5 
0.2 



"IQ 



0.18 
0.08 
0.026 
0.015 



£r 



10.0 
12.8 

14.5 
18.0 



TABLE IV: Heavy and light hole masses long with 
background dielectric constants for various III-V 
semiconductors — '. 
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FIG. 8: The imaginary part Im(e-"^'*(q,c^)) of the the RPA 
dielectric function as a function of wave number q and energy 
Tujj for various semiconductor systems at a total hole density 
ofn = 0.35nm-^ 
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FIG. 9: The modulus |e^^'*(g,aj)| of the RPA dielectric 
function as a function of wave number q and energy huj 
for various semiconductor systems at a total hole density of 
n = 0.35nm~^. The dark areas indicate zeros of the dielectric 
function corresponding to plasmon excitations. 
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FIG. 10: The real (top panels) and imaginary (bottom pan- 
els) part of the free polarizability xo(9*, i^) as a function of 
frequency at different wave vectors for the same semiconduc- 
tor systems as before (cf. table |IV|. 



simplify considerably tcP^ll 
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2kh 
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2tt-^K^^ \2kh 



-L 



kh + h 
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kh + ki 
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2Tr^h^'^^ \2k, 



kh + h 
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2Tr^h? "" \kh + ki^ 
where L(x) is the so-called Lindhard correction, 






1 



1-x 



(31) 



(32) 



and the function H is defined as 



H{x) 



l/x ^ 

dy— In 



y 

-vOO 



1 + 2/ 



1-2/ 



II _ V°° ^''"^' Irl < 1 
Z_/n=0 (2n+l)2 FI - 1 



(33) 



Details of the derivation of the above result can be found 
in appendix I A 1 a[ Note that the static polarizability can 
entirely be expressed in terms of the arguments k/2kh, 
k/2ki, and k/kh + ki with the latter one being the har- 
monic mean of the two former. In the case nih — mi 




FIG. 11: Left panel: The static free polarizability Xoi<l, 0) for 
the same choice of parameters as in Figs. |2)6[ Right panels: 
Xo{<l, 0) and e^^ {q, 0) the same III-V semiconductor systems 
as in Figs. [Tflol (cf. table |lvf. 



(i.e. kh = ki =: kp) one obtains the usual result 
Xo(9jO) = —D{eF)L{q/2kp) for charge carriers without 
spin-orbit coupling where D{e) is the density of statea^^. 
For mil y^ mi, however, the static polarizability ( |31[ ) has 
a clearly more complicated structure. Fig. ITT] displays 
the static free polarizability and dielectric function for 
the systems discussed above. In particular, the data in 
the left panel at fixed mh+mi — mp shows that the static 
polarizability develops richer features with increasing dif- 
ference in heavy and light hole mass. 

Moreover, in the long- wave approximation xo (?: 0) ~ 
Xo(0,0) one recovers the usual Thomas-Fermi (TF) 
screening, 



Hg,0)«l-4^/g2 



(34) 



Qtf 



(eV£r£o)3n/(2e/). 

As discussed in Ref.'^, the full screened potential of a 
pointlike probe charge Q, 



$(r) 



(27r)3 



'3„ ErEoq^ iqf 



d-'q 



e 



RPA 



(q) 



(35) 



can conveniently be approximated using Lighthill's the- 
orenpSl as 

Hr) « "^0oo(2A:^, r) + '^(/.^(2fc;, r) (36) 

Too mo 



where 



^ociq-.r) 



Q 



cos{qr) 



47r£oao7r(£^eflPA(^))2 (^^)3 



(37) 



and uq = ATTS^h? / {niQe^) being the usual Bohr radius. 
As a result, a beating of Friedel oscillations between the 



two wave numbers 2fc^/; (but not k = [kh + ki)/2) takes 
place^. This beating is a peculiarity of the holes residing 
in the p-type valence band and should be observable via 
similar scanning tunneling microscopy techniques as used 
in metals^^ and n-doped semiconductors^^. Moreover, as 
theoretical studies have revealed, such oscillations can 
have a profound impact on the magnetic properties of 
ferromagnetic semiconductors by giving way to the pos- 
sibility of noncoUinear magnetic orderin g^^-^^' . 



B. Limit of large frequencies 

In the regime of large frequencies and small wave vec- 
tors, one can expand the denominators in Eq. ( |10p as- 
suming huj » eh/i{q) and fno » {Hkii/i/m,fi^i)hq. The 
result within the two leading orders reada^ESl 



.RPA 



(q^^) 
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2;i2 1 / 1 1 
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mh 
1 
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9^ (fcf. 






Y2^' (4 + kf) 



ijJ^ ErEQTT^ 



21 



200^ 
3 2 



1 
56 

1 

1 
rrih 



1 

rrih 

1 

m? 



1 

mi 



{kl 



(kl kj) 



kf) 



1 

mi 



1,5 



kf 



(38) 



For mil = mi the first three lines of the above expression 
reproduce again the standard textbook result^ while all 
other terms vanish in this limit. On the other hand, if 
TO./,. ^ mi, contributions in order l/w^ occur that are in- 
dependent of the wave vector q. Such terms are absent 
in the case of the standard electron gas where the contri- 
butions of order l/w^" are at least of order g^"~^ in the 
wave vector^. The technical reason why such contribu- 
tions are present for the hole gas is that the expression 
E\2{k + q) — £\i{k) in Eq.(lO) contains for |Ai| ^ IA2I an 



additive term which is independent of k (and vanishes for 
mil = rni). As a consequence, although the result (38) 



is the valid high frequency expansion of the dielectric 
function, it is not possible to obtain from it a reliable 
expression for the plasmon dispersion uj{q) defined by 
E^^^'^{q,(jj{q)) = 0. This statement holds even for the 
long- wavelength plasma frequency uj{q = 0) and is due 
to the fact that in any order in (1/w^) the prefactor in 
the expansion contains contributions being of low order 
(including zeroth order) in q. As an example, relying on 
the expansion (38) being of up to quartic order in l/cj, 
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the condition £^^^{q,uj{q)) = translates tcP^SI 

Here we have defined 

«0'-— !^fz^ + Z^V (40) 



£^£0 2 \mh mi 



and the density-independent coefScients u, w, and w are 
given by 



-Q{mh,mi) 



(37r2)i/3 (• 



1 1 

14 V m/i TO/ 



3/2 , 3/2 



2/3 



)(„.;'' -mP). (41) 



(4 + ^)H"""'n. (^^' 



2Q(to/i,TO/) / 1 1 \ / 5/2 , 5/2 



w 



3Q{mh,mi) 
107r2 



( 




5/2 
TO^' 

5/2 \ 



5/2 



mt 



with the common prefactor 



Q{mh,mi) = 



t^ruQ 



(3^2) 



2^5/3 



1 



4) (^ 



3/2 



3/2 



5/3 



(43) 



(44) 



Clearly, the coefficients u and w vanish for mh = mi while 
from V one recovers usual textbook result for an electron 
gas without spin-orbit coupling^!. However, if rrih and 
TO,; differ substantially, the neglected contributions oc- 
curring in higher order in the inverse frequency but be- 
ing independent of or of low order in the wave vector can 
substantially modify the plasmon excitations. This can 
even affect the plasma frequency w(g = 0) at zero wave 
vector: E.g. for the parameters of GaAs (cf. table IV) 
and a total hole density of n 
from Eq. ^ huj{q = 0) 



0.35nm~'^ one obtains 
0.8eV, in contrast to the 
value of slightly less than 0.3eV found from a full evalua- 
tion of the dielectric function (cf. Fig. l9|top right panel) 
which is also in accordance with Ref. . In summary, 
although the expansion [38] is the correct description of 
the dielectric function in the limit of large frequencies, 
it does not lead to reliable expressions for the plasmon 
dispersion. This is due to peculiarities of the expansion 
occurring for rrih ^ TO;. On the other hand, as seen from 
e.g. Fig.|9j the interplay between the hole mass difference 



and the background dielectric constant leads to plasmon 
excitations patterns which differ dramatically from the 
textbook case of the jellium model. This alternations, 
however, are not accurately described by expressions of 
the type ( 39 1 , in contrast to an earlier approach where 



results for the full wave vector and frequency dependence 
of the dielectric function were not available yetP. 



V. CONCLUSIONS AND OUTLOOK 

We have investigated the RPA dielectric function of the 
homogeneous semiconductor hole liquid in p-doped bulk 
III-V zinc-blende semiconductors. The single-particle 
physics of the hole system is modeled by Luttinger's 
four-band Hamiltonian in its spherical approximation. 
Regarding the Coulomb-interacting hole liquid, the full 
dependence of the zero-temperature dielectric function 
on wave vector and frequency has been explored. The 
imaginary part of the dielectric function is analytically 
obtained in terms of complicated but fully elementary 
expressions, while in the result for the real part nonele- 
mentary one-dimensional integrations remain to be per- 
formed. The correctness of these two independent calcu- 
lations is checked via Kramers-Kronig relations. 

The mass difference between heavy and light holes, 
along with variations in the background dielectric con- 
stant, leads to dramatic alternations in the plasmon ex- 
citation pattern, and generically, two plasmon branches 
can be identified. These findings are the result of the 
evaluation of the full dielectric function and are not ac- 
cessible via a high-frequency expansion. In the static 
limit a beating of Friedel oscillations between the Fermi 
wave numbers of heavy and light holes occurs. 

Regarding future developments, possible extensions of 
the present work could include the implementation of 
more general single-particle Hamiltonians modeling the 
band structure. For instance, one could drop the spheri- 
cal approximation to the Luttinger Hamiltonian and con- 
sider parameters 72 7^ 73. However, such a reduction 
of the full rotational invariance to tetrahedral symmetry 
might prohibit analytical progress as achieved here. How- 
ever, we do not expect drastic effects from such a gener- 
alization, in particular not since for thegeneric material 
GaAs 72, 73 are very close to each othepSS. Moreover, our 
results obtained here for the spherically symmetric 4x4 
valence band Hamiltonian agree, where overlapping, very 
reasonably with findings in Ref.l^ where a more compli- 
cated 8x8 band structure model was evaluated numeri- 
cally. 

Having in mind ferromagnetic semiconductors such as 
Mn-doped GaAs, another obvious extension would be a 
coupling to the hole spins by a homogeneous Zeeman- 
type field mimicking the magnetization of the Mn ions. 
A technical difficulty here lies in the fact that the result- 
ing single-particle Hamiltonian cannot be diagonalized 
any more in a convenient analytical fashion. However, 
analytical progress might still be possible if the Zeeman 
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coupling is treated as a perturbation. This strategy leads 
of course to also consider the spin susceptibility, i.e. spin- 
spin response function. For a standard jellium systems of 
spin-1/2 fermions without spin-orbit coupling, this quan- 
tity is, up to constant prefactors, identical to the free 
electrical polarization^. For the 4x4 hole system studied 
here, however, this simple relation is rendered invalid by 
the larger spin length and, more importantly, the pres- 
ence of manifest spin-orbit coupling. Thus, a study of 
the spin susceptibility in a similarly analytical fashion as 
done here for the electric polarizability and the dielectric 
function appears also to be useful. 

Finally, from the point of view of general many-particle 



physics, the inclusion of so-called local many-body field 
factors would be an important step towards correlations 
beyond RPAQJ. The practical treatment of such local field 
factors in the presence of strong spin-orbit coupling, how- 
ever, is still in its infancy. 
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Appendix A: Calculation of the free polarizability 
1. The real part 

With the help of the Dirac identity 

— — - =V- -iTi5{x) 

X + l() X 

the contributions to the real part of the free polarizability of the hole gas read 



(Al) 



Rhh{q,i^) = 



Rhi{q,ijj) = 
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dx\ 6- 3—— — — 
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1 + {2k/q)x + (1 - nii/m,h)k'^/q'^ - 2miuj/{hq'^' 
The integration over the polar variable x can be performed by applying the identity 
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(A4) 



for a = —2mhUj/(hq^) and a = (1 — mi/mh)k'^ /q'^ — 2miuj/{hq^). Adding both contributions, and introducing a 
dimensionless radial integration variable y — 2k/ q, the result can be formulated as 
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Now, by rearranging the terms in the integrand and performing elementary integrations, one obtains the result ( 14 1 
We note that also the integral in the third line of Eq. ( 14 1 is elementary but lengthy. 
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whereas all other integrals in Eq. ( 14 1 cannot be expressed via elementary functions. 
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a. Static limit 



In the static limit w — > one obtains from Eq. (14 1 (or, alternatively, Eq. (A5)) 
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(A6) 



where we have split a part of the logarithmic terms in the integrand by introducing y± := 2 / {\± ^J mi / mh) ■ In order to 
simplify the above expression we first note that the first term of the r.h.s can be written as —mhL{q/{2kh))/{T^fif' using 
the Lindhard correction (32). Regarding the terms in the second and third line involving a summation over /i = ±, 



the interchange h O / leads to y± o y± :— 2/(1 ± yjmhjmi) fulfilling 2ki/[qy±) = ±2kh/{qy±) and mhy± — raiy\. 
From these observations it is easy to see that the terms with /i = — cancel against corresponding expressions in 
Rii{q, 0) + Rih{q, 0), and only the terms with /i = + (being invariant under h -{^ I) contribute to Xo(9, 0). The first of 
these contributions (second line in Eq. (A6)) can be expressed again via the Lindhard correction, while the integrals 
in the third line involve the function H[x) defined in Eq. (33). Finally, the last line of Eq. (A6) can be evaluated as 
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(A7) 



where the last term on the r.h.s. is odd under h -(->■ I and cancels agianst an analogous contribution in Rii(q,0) 
Rih{q,0). Now summing all expressions one ends up with the result (31) for the free polarizability Xo{q,0)- 



2. The imaginary part 

a. Ihh{q,^) and Iu{q,uj) 



Using the Dirac identity (Al) and performing the angular integrations, Ihhiq^^) can be expressed as 
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(A8) 
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where Q{x) denotes the Heaviside step function. Ob- 
viuosly, the step functions occurring in the above ex- 
pression define the lower integration bound, and the 
pertaining discussion parallels the arguments appropri- 
ate for the standard textbook case of a spinless Jellium 
model^ ''^. However, for the sake of completeness, and 
in order to point out important differences regarding the 
remaining quantities Ihi{q,(^) and Iih{q,uj) to be ana- 
lyzed below, let us briefly go into some details. Since 
//i^((f, — cj) = — //i/i((f, w) it is sufficient to concentrate on 
u! > 0. Then, if the first step function in (A8) leads to 



a non-zero contribution (i.e. has a positive argument for 
some k G [0, kh]), this holds also for the second step func- 
tion. Thus, a necessary and sufhcient condition for both 
step function to contribute is kh — \q/2 + mhUo/{hq)\ > 0, 
which is equivalent to 



hq Ehiq) > huj 

mil 



(A9) 



and can for non-negative frequencies only be fulfilled for 
q < 2kh. The last two inequalities define region I in 
table|l] and the corresponding expression ( 15 1 is obtained 
by elementary integration. 

Let us now turn to the case where only the second 
step function con tribu tes, i.e. kh — \q/2 — mhOJ / {hq)\ > 



while inequality (A9) is violated. Assuming hw > eh{q) 
we arrive at the condition 

hq^^ + Ehiq) >hLO> hq^^ - eh{q) , (AlO) 

ruh rrih 



while the opposite assumption, huj < Sh{q), leads to 

(All) 



Hkh , , 

nuj > —hq h Sh(q) ■ 

rrih 



The latter inequality is only a nontrivial condition if its 
r.h.s. is non- negative which is equivalent to g > 2kh- 



In summary, the second step function in expression ( A8 ) 



contributes while the first one yields zero if, and only 
if, (i) q ^ 2fc;i and inequality (AlO) is fulfilled, or (ii) 
q > 2kh and 



,- frkh , , . 

nq V Sh\q) > nuj > 

ruh 



-hq- 



hkh 
rrih 



+ eh{q). (A12) 



The above conditions define region II in table IT] and and 
the corresponding expression (16) follows again from el- 
ementary integration. 

The remaining quantity Iii{q,uj) is obtained from the 
above results via the replacement h i~^ I. 



b. Ihi{q,uj) 

After performing the angular integrations, Ihi{q, i-^) can 
be formulated in the form (17) with 
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(A13) 



We now have to discuss under which circumstances the 
step functions in the above expression lead to nonzero 
contributions. This is more complicated than in the pre- 
vious case since the arguments of the step functions de- 
pend quadratically (and not only linearly) on the inte- 
gration variable k. The condition 



where we have defined 



e /c 



is equivalent to 



mi \ k . miuj 
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(A15) 
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(A16) 



(A17) 



Here and in what follows the upper (lower) sign refers 
always to I^i{q,uj) {I^i{q,uj)). Note that the step func- 
tion in I'^i{q,uj) can, for non-negative frequencies, only 
be nonzero if this is also the case for I^i{q,uj). 

Since m; < ruh we clearly have a > 0, and the second 
inequality in (A15) requires b± > which leads for the 
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upper sign to the condition 



huj < 



1 



-£h{q)^ 



(A18) 



1 - mi/mh 
Moreover, an elementary discussion of the inequahties 



(A15) yields the following lower and upper boundaries 



for the integral (A13) after resolving the step function, 

l± = min{|a- v/fc^|,fc,,} , (A19) 

/± = mini a+ ^/b±,kh\ , (A20) 

with < i± < ^± < fc/i- Nonzero contributions occur 
only for l^. < kfi. For the upper sign, the condition l^. > 
kh implies 



hq 



hkh 



^liq)- 



rrih 



— 1 \ Ef < huj 



(A21) 



mi \mi 

which is, for non-negative w, a nontrivial statement on if 

(l - \Jrailmi^ K < q < {^ + \Jmilmi}^ ky, . (A22) 



Conversely, inequality (A21) implies l_j^ > kh only for 



a < kh, while in the case a > kh, i.e. 
q < {I ~mi/mh)kh.. 



(A23) 



it follows kh < a + •v/fe+, and the inequality (A18) limits 



the re gion o f non zero I'f^i{q,uj). The inequalities (A23l, 
( |A21[ ), and ( |A18| ) define region II in table |ll] with the 
integration bounds being l_j^ — \a — \/b 



k^ and 
b+ =: kh as defined in Eqs. (|21l),(|22l). On the 



-r+ 



1+ = a 

other hand, inequality (|A22|) along witE~the"nlgation of 



(A21) 



hq^ - sM) -{—-l]sf>nu:, (A24) 

mi \mi J ' 

define region I in table II Here again /_,_ = Kh, and in- 
equality (A24| further implies /+ = kh- Note that the 



condition (A18) is always fulfilled if (A21) is valid since 

£h{q) 



hq 



hkh 



mi 
^ 0< 



£iiq)- 
h^ 

2mi 



mi 
mi 



\^—\kh-c 



mi/nih 



mi 
mh 



Moreover, the upper and lower boundary of region II in- 
tersect each other in the g-w-plane ai q = (1 — mi/mh)kh 
with identical tangent. Finally, the corresponding contri- 
butions to I'^i {q, w) in r egions I and II are obtained by ele- 



mentary integration of (A13) and given in Eqs. (|19[),(20|. 
Let us now turn to the lower sign case /^;((f,tj). The 
condition /_ < kh implies for 6_ > a^ (<=> hu > ei{q)) 



hq- 



hkh 
mi 



^i{q) 



mh 
nil 



i\ Ef >hijj . 



(A25) 



In the opposite case 5_ < a? (<^ hu < £/((?)) the inequal- 
ity I- < kh does not lead to any further restriction on 
the frequency for kh > a (<=> (1 — mi/mh)kh > q), while 
for kh < a one finds 



hw > —hq h ei{q) 

mi 



mh 
mi 



1 



(A26) 



The latter statement is a nontrivial requirement only if 



q> [1 + \/mi/mh] kh 



(A27) 



which also ensures (1 — mi/mh)kh < q {-^ kh < a). The 
inequalities (A25l and (A26) are necessary and sufficient 



conditions for If^i{q,uj) in (Al3l to be nonzero. In this 
case the lower integration bound is {_ = \a— \/hZ\ =: k^^ 



and again given explicitly in Eq. (21) 



Moreover, straightforward inspection shows that the 
upper integration boundary is I- = kh provided inequal- 



ity (A26l ( but not necessarily (A27)) is fulfilled, while 
otherwise (requiring q < (1 — \J fn^i / m,h)kh) we have 



6_ =: fc,j as given in Eq. 




The corre- 
([24| follow 



l_ = a - 

spending contributions G_(. . . ) in Eqn 
again from elementary integration. The different regions 
in the q-w-plane are summarized in tabic |ll] and illus- 
trated in Fig. [T] 



c. Iih{q,i^) 

The contr ibution Iih{q,u!) can formally be expressed 
by Eq. ( A13l via the replacement h ^ I. Thus one needs 



to discuss the condition 
mh 



e k- 



or, equivalently. 



1- 



mi 



± 



nihU} 



2q hq 



= 1, 



{k- 



> di 



A (fc - cf < d± 



with 



d± 



c' + 



mh/mi - 1 ' 

rrih/mi - 
(mh/mi)q^ 
{mh/mi - 1)2 



q nihuj 

2 hq 

2mhi-tj/h 

mh/mi - 1 



(A28) 



(A29) 



(A30) 



(A31) 



Note that, differently from the previous cases, Ii^{q,uj) 
is not necessarily zero if Il^(q,ui) vanishes, because 1 — 
nfih/mi < 0. On the other hand, this inequality ensures 
c > 0, and from the condition d-t > we find for the 
lower case 



htj < 



1 



mh/mi ~ 1 



ei{q)- 



(A32) 
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Similarly to the previous case, the inequalities (A29) and 



lead to the following lower and upper boundaries within 
the interval [0, fc;], 



l± 


= min i |c — 


Vd±\,ki\ , 


(A33) 


l± 


= min ■ c + y 


/d^,h} , 


(A34) 



with nonzero contributions being possible only for l^. < 
ki. For the upper sign case this condition is equivalent 
to 

hq^ - suiq) -(—-l]ef> hi, (A35) 

which can only be fulfilled if 

q<(l + VmJ^i) h . (A36) 

Thus, in the above case, the lower integrati on b ound is 



d+l =■ k^ and given explicitly in (27 1. More- 
over, an again straightforward discussion shows that the 



upper integration bound is given by /+ 
(cf. Eq. ^) provided 



C+ -v/Zi 



hw < 



. hki f mi 
-hq eh{q) 1 £/ > 



which is only possible for 

q < (-1 + y^mh/mA ki . 



(A37) 



(A38) 



The above inequalities ( A37), (A38[) define region II in ta- 



ble 



: III and Fig. fl] while region I is defined by ( A35 1 , ( A36 ) 
and the negation of (A37). Here the upper integration 
bound is 1+ — ki. The corresponding express ion s for 



Ii^{q,uj) in regions I, II are given in Eqs. (25), (26), re- 
spectively. 

Regarding Ij^{q,uj), considerations analogous to the 
ones for /^^((f, oj) show that the condition /_ < ki is for 
c > ki (<^ q > (mh/mi — l)fc;) equivalent to 



hq \-ehiq) 

rrih 



■mi 
mfi 



1 \ Ef > hcu 



(A39) 



Hco > -hq^ + Shiq) + f ^ - 1 ) £/ , (A40) 



where latter inequality poses a nontrivial requirement 
only if 



q>[l + y/nih/mA kh 



(A41) 



The inequahties ( |A39[ ), ( |A40[ ), and \Ml\ define region 
III in table III Here the lower integration bound is /_ = 



|c— ■\/(i-| =: fc; (cf. Eq. (27)), and the upper integration 
bound turns out to be always I_ = /c/. 
For c < ki, i.e. 



q < {rrih/mi - 1) /c/ 



(A42) 



however, there is, similar to the case of /^^((f, w), another 
way of fulfilling the condition /„ < ki. The corresponding 
region IV is defined by the inequalities ( A32 1, ( A42 ), and 



hLu>hq—+eh{q)+(—-l]sf, (A43) 



and the integration bounds are l^^k^ , ?_ = 



/dZ=: 



^l : 



We note that the inequality (A39) generally implies 
the fulfillment of ( A32 ) since 



hq \-eh{q) 

■m-h 



<^ 



0< 



mi ^ 
mh 
f mh 
2mh V V "^; 



1 £/< 



£i{q) 



mh/mi - 1 
-l]ki-q' 



Moreover, similarly as in the case of I'^i{q,uj), the upper 
and the lower boundary of region IV intersect each other 
at q = (rnh/m,i — l)ki with identical tangent. 
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